pH-dependent reaction triggering in PmHMGR crystals for time-resolved crystallography.

Serial crystallography and time-resolved data collection can readily be employed to investigate the catalytic mechanism of Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl (HMG)-coenzyme-A (CoA) reductase (PmHMGR) by changing the environmental conditions in the crystal and so manipulating the reaction rate. This enzyme uses a complex mechanism to convert mevalonate to HMG-CoA using the co-substrate CoA and cofactor NAD+. The multi-step reaction mechanism involves an exchange of bound NAD+ and large conformational changes by a 50-residue subdomain. The enzymatic reaction can be run in both forward and reverse directions in solution and is catalytically active in the crystal for multiple reaction steps. Initially, the enzyme was found to be inactive in the crystal starting with bound mevalonate, CoA, and NAD+. To observe the reaction from this direction, we examined the effects of crystallization buffer constituents and pH on enzyme turnover, discovering a strong inhibition in the crystallization buffer and a controllable increase in enzyme turnover as a function of pH. The inhibition is dependent on ionic concentration of the crystallization precipitant ammonium sulfate but independent of its ionic composition. Crystallographic studies show that the observed inhibition only affects the oxidation of mevalonate but not the subsequent reactions of the intermediate mevaldehyde. Calculations of the pKa values for the enzyme active site residues suggest that the effect of pH on turnover is due to the changing protonation state of His381. We have now exploited the changes in ionic inhibition in combination with the pH-dependent increase in turnover as a novel approach for triggering the PmHMGR reaction in crystals and capturing information about its intermediate states along the reaction pathway.

Heterogeneity in Lowe Syndrome: Mutations Affecting the Phosphatase Domain of OCRL1 Differ in Impact on Enzymatic Activity and Severity of Cellular Phenotypes.

Lowe Syndrome (LS) is a condition due to mutations in the OCRL1 gene, characterized by congenital cataracts, intellectual disability, and kidney malfunction. Unfortunately, patients succumb to renal failure after adolescence. This study is centered in investigating the biochemical and phenotypic impact of patient's OCRL1 variants (OCRL1VAR). Specifically, we tested the hypothesis that some OCRL1VAR are stabilized in a non-functional conformation by focusing on missense mutations affecting the phosphatase domain, but not changing residues involved in binding/catalysis. The pathogenic and conformational characteristics of the selected variants were evaluated in silico and our results revealed some OCRL1VAR to be benign, while others are pathogenic. Then we proceeded to monitor the enzymatic activity and function in kidney cells of the different OCRL1VAR. Based on their enzymatic activity and presence/absence of phenotypes, the variants segregated into two categories that also correlated with the severity of the condition they induce. Overall, these two groups mapped to opposite sides of the phosphatase domain. In summary, our findings highlight that not every mutation affecting the catalytic domain impairs OCRL1's enzymatic activity. Importantly, data support the inactive-conformation hypothesis. Finally, our results contribute to establishing the molecular and structural basis for the observed heterogeneity in severity/symptomatology displayed by patients.

Computational Study of Base-Catalyzed Thiohemiacetal Decomposition in Pseudomonas mevalonii HMG-CoA Reductase.

Thiohemiacetals are key intermediates in the active sites of many enzymes catalyzing a variety of reactions. In the case of Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl coenzyme A reductase (PmHMGR), this intermediate connects the two hydride transfer steps where a thiohemiacetal is the product of the first hydride transfer and its breakdown forms the substrate of the second one, serving as the intermediate during cofactor exchange. Despite the many examples of thiohemiacetals in a variety of enzymatic reactions, there are few studies that detail their reactivity. Here, we present computational studies on the decomposition of the thiohemiacetal intermediate in PmHMGR using both QM-cluster and QM/MM models. This reaction mechanism involves a proton transfer from the substrate hydroxyl to an anionic Glu83 followed by a C-S bond elongation stabilized by a cationic His381. The reaction provides insight into the varying roles of the residues in the active site that favor this multistep mechanism.

Targeting Enterococcus faecalis HMG-CoA reductase with a non-statin inhibitor.

HMG-CoA reductase (HMGR), a rate-limiting enzyme of the mevalonate pathway in Gram-positive pathogenic bacteria, is an attractive target for development of novel antibiotics. In this study, we report the crystal structures of HMGR from Enterococcus faecalis (efHMGR) in the apo and liganded forms, highlighting several unique features of this enzyme. Statins, which inhibit the human enzyme with nanomolar affinity, perform poorly against the bacterial HMGR homologs. We also report a potent competitive inhibitor (Chembridge2 ID 7828315 or compound 315) of the efHMGR enzyme identified by a high-throughput, in-vitro screening. The X-ray crystal structure of efHMGR in complex with 315 was determined to 1.27 Å resolution revealing that the inhibitor occupies the mevalonate-binding site and interacts with several key active site residues conserved among bacterial homologs. Importantly, 315 does not inhibit the human HMGR. Our identification of a selective, non-statin inhibitor of bacterial HMG-CoA reductases will be instrumental in lead optimization and development of novel antibacterial drug candidates.

Essential amino acids in the Plant-Conserved and Class-Specific Regions of cellulose synthases.

The Plant-Conserved Region (P-CR) and the Class-Specific Region (CSR) are two plant-unique sequences in the catalytic core of cellulose synthases (CESAs) for which specific functions have not been established. Here, we used site-directed mutagenesis to replace amino acids and motifs within these sequences predicted to be essential for assembly and function of CESAs. We developed an in vivo method to determine the ability of mutated CesA1 transgenes to complement an Arabidopsis (Arabidopsis thaliana) temperature-sensitive root-swelling1 (rsw1) mutant. Replacement of a Cys residue in the CSR, which blocks dimerization in vitro, rendered the AtCesA1 transgene unable to complement the rsw1 mutation. Examination of the CSR sequences from 33 diverse angiosperm species showed domains of high-sequence conservation in a class-specific manner but with variation in the degrees of disorder, indicating a nonredundant role of the CSR structures in different CESA isoform classes. The Cys residue essential for dimerization was not always located in domains of intrinsic disorder. Expression of AtCesA1 transgene constructs, in which Pro417 and Arg453 were substituted for Ala or Lys in the coiled-coil of the P-CR, were also unable to complement the rsw1 mutation. Despite an expected role for Arg457 in trimerization of CESA proteins, AtCesA1 transgenes with Arg457Ala mutations were able to fully restore the wild-type phenotype in rsw1. Our data support that Cys662 within the CSR and Pro417 and Arg453 within the P-CR of Arabidopsis CESA1 are essential residues for functional synthase complex formation, but our data do not support a specific role for Arg457 in trimerization in native CESA complexes.

Microsecond timescale MD simulations at the transition state of PmHMGR predict remote allosteric residues.

Understanding the mechanisms of enzymatic catalysis requires a detailed understanding of the complex interplay of structure and dynamics of large systems that is a challenge for both experimental and computational approaches. More importantly, the computational demands of QM/MM simulations mean that the dynamics of the reaction can only be considered on a timescale of nanoseconds even though the conformational changes needed to reach the catalytically active state happen on a much slower timescale. Here we demonstrate an alternative approach that uses transition state force fields (TSFFs) derived by the quantum-guided molecular mechanics (Q2MM) method that provides a consistent treatment of the entire system at the classical molecular mechanics level and allows simulations at the microsecond timescale. Application of this approach to the second hydride transfer transition state of HMG-CoA reductase from Pseudomonas mevalonii (PmHMGR) identified three remote residues, R396, E399 and L407, (15-27 Å away from the active site) that have a remote dynamic effect on enzyme activity. The predictions were subsequently validated experimentally via site-directed mutagenesis. These results show that microsecond timescale MD simulations of transition states are possible and can predict rather than just rationalize remote allosteric residues.

Visualizing the enzyme mechanism of mevalonate diphosphate decarboxylase.

Mevalonate diphosphate decarboxylases (MDDs) catalyze the ATP-dependent-Mg2+-decarboxylation of mevalonate-5-diphosphate (MVAPP) to produce isopentenyl diphosphate (IPP), which is essential in both eukaryotes and prokaryotes for polyisoprenoid synthesis. The substrates, MVAPP and ATP, have been shown to bind sequentially to MDD. Here we report crystals in which the enzyme remains active, allowing the visualization of conformational changes in Enterococcus faecalis MDD that describe sequential steps in an induced fit enzymatic reaction. Initial binding of MVAPP modulates the ATP binding pocket with a large loop movement. Upon ATP binding, a phosphate binding loop bends over the active site to recognize ATP and bring the molecules to their catalytically favored configuration. Positioned substrates then can chelate two Mg2+ ions for the two steps of the reaction. Closure of the active site entrance brings a conserved lysine to trigger dissociative phosphoryl transfer of γ-phosphate from ATP to MVAPP, followed by the production of IPP.

Conformational Dynamics in the Binding-Protein-Independent Mutant of the Escherichia coli Maltose Transporter, MalG511, and Its Interaction with Maltose Binding Protein.

MalG511 is a genetically selected binding-protein-independent mutant of the Escherichia coli maltose transporter MalFGK2, which retains specificity for maltose and shows a high basal ATPase activity in the absence of maltose binding protein (MBP). It shows an intriguing biphasic behavior in maltose transport assays in the presence of MBP, with low levels of MBP stimulating the activity and higher levels (>50 µM) inhibiting the transport activity. Remarkably, the rescuing effect of the MBP suppressor mutant, MBPG13D, turns it into an attractive model for studying regulatory mechanisms in the ABC transporter superfamily. It is hypothesized that the special characteristics of MalG511 result from mutations that shift its equilibrium toward the transition state of MalFGK2. We tested this hypothesis by using site-directed spin labeling in combination with electron paramagnetic resonance spectroscopy, which showed conformational changes in MalG511 and its interaction with MBP and MBPG13D during its catalytic cycle. We found that MalG511 utilizes the same alternate access mechanism as MalFGK2, including all three open, semi-open, and closed states of the MalK dimer, to transport maltose across the membrane. However, the equilibrium of this mutant is shifted toward the semi-open state in its resting state and interacts with MBP with high affinity, providing an explanation for the inhibition of MalG511 by MBP at higher concentrations. In contrast, the mutant binding protein, MBPG13D, interacts with lower affinity and could restore MalG511 to a normal catalytic cycle.

Mevalonate 5-diphosphate mediates ATP binding to the mevalonate diphosphate decarboxylase from the bacterial pathogen Enterococcus faecalis.

The mevalonate pathway produces isopentenyl diphosphate (IPP), a building block for polyisoprenoid synthesis, and is a crucial pathway for growth of the human bacterial pathogen Enterococcus faecalis The final enzyme in this pathway, mevalonate diphosphate decarboxylase (MDD), acts on mevalonate diphosphate (MVAPP) to produce IPP while consuming ATP. This essential enzyme has been suggested as a therapeutic target for the treatment of drug-resistant bacterial infections. Here, we report functional and structural studies on the mevalonate diphosphate decarboxylase from E. faecalis (MDDEF). The MDDEF crystal structure in complex with ATP (MDDEF-ATP) revealed that the phosphate-binding loop (amino acids 97-105) is not involved in ATP binding and that the phosphate tail of ATP in this structure is in an outward-facing position pointing away from the active site. This suggested that binding of MDDEF to MVAPP is necessary to guide ATP into a catalytically favorable position. Enzymology experiments show that the MDDEF performs a sequential ordered bi-substrate reaction with MVAPP as the first substrate, consistent with the isothermal titration calorimetry (ITC) experiments. On the basis of ITC results, we propose that this initial prerequisite binding of MVAPP enhances ATP binding. In summary, our findings reveal a substrate-induced substrate-binding event that occurs during the MDDEF-catalyzed reaction. The disengagement of the phosphate-binding loop concomitant with the alternative ATP-binding configuration may provide the structural basis for antimicrobial design against these pathogenic enterococci.

The nuclear mitotic apparatus protein NuMA controls rDNA transcription and mediates the nucleolar stress response in a p53-independent manner.

The nuclear mitotic apparatus protein, NuMA, is involved in major cellular events such as DNA damage response, apoptosis and p53-mediated growth-arrest, all of which are under the control of the nucleolus upon stress. Proteomic investigation has identified NuMA among hundreds of nucleolar proteins. Yet, the precise link between NuMA and nucleolar function remains undetermined. We confirm that NuMA is present in the nucleolus and reveal redistribution of NuMA upon actinomycin D or doxorubicin-induced nucleolar stress. NuMA coimmunoprecipitates with RNA polymerase I, with ribosomal proteins RPL26 and RPL24, and with components of B-WICH, an ATP-dependent chromatin remodeling complex associated with rDNA transcription. NuMA also binds to 18S and 28S rRNAs and localizes to rDNA promoter regions. Downregulation of NuMA expression triggers nucleolar stress, as shown by decreased nascent pre-rRNA synthesis, fibrillarin perinucleolar cap formation and upregulation of p27kip1, but not p53. Physiologically relevant nucleolar stress induction with reactive oxygen species reaffirms a p53-independent p27kip1 response pathway and leads to nascent pre-rRNA reduction. It also promotes the decrease in the amount of NuMA. This previously uncharacterized function of NuMA in rDNA transcription and p53-independent nucleolar stress response supports a central role for this nuclear structural protein in cellular homeostasis.

Rice Cellulose SynthaseA8 Plant-Conserved Region Is a Coiled-Coil at the Catalytic Core Entrance.

The crystallographic structure of a rice (Oryza sativa) cellulose synthase, OsCesA8, plant-conserved region (P-CR), one of two unique domains in the catalytic domain of plant CesAs, was solved to 2.4 Å resolution. Two antiparallel α-helices form a coiled-coil domain linked by a large extended connector loop containing a conserved trio of aromatic residues. The P-CR structure was fit into a molecular envelope for the P-CR domain derived from small-angle X-ray scattering data. The P-CR structure and molecular envelope, combined with a homology-based chain trace of the CesA8 catalytic core, were modeled into a previously determined CesA8 small-angle X-ray scattering molecular envelope to produce a detailed topological model of the CesA8 catalytic domain. The predicted position for the P-CR domain from the molecular docking models places the P-CR connector loop into a hydrophobic pocket of the catalytic core, with the coiled-coil aligned near the entrance of the substrate UDP-glucose into the active site. In this configuration, the P-CR coiled-coil alone is unlikely to regulate substrate access to the active site, but it could interact with other domains of CesA, accessory proteins, or other CesA catalytic domains to control substrate delivery.

Biochemical characterization of essential cell division proteins FtsX and FtsE that mediate peptidoglycan hydrolysis by PcsB in Streptococcus pneumoniae.

The FtsEX:PcsB complex forms a molecular machine that carries out peptidoglycan (PG) hydrolysis during normal cell division of the major respiratory pathogenic bacterium, Streptococcus pneumoniae (pneumococcus). FtsX is an integral membrane protein and FtsE is a cytoplasmic ATPase that together structurally resemble ABC transporters. Instead of transport, FtsEX transduces signals from the cell division apparatus to stimulate PG hydrolysis by PcsB, which interacts with extracellular domains of FtsX. Structural studies of PcsB and one extracellular domain of FtsX have recently appeared, but little is known about the biochemical properties of the FtsE ATPase or the intact FtsX transducer protein. We report here purifications and characterizations of tagged FtsX and FtsE proteins. Pneumococcal FtsX-GFP-His and FtsX-His could be overexpressed in Escherichia coli without toxicity, and FtsE-His remained soluble during purification. FtsX-His dimerizes in detergent micelles and when reconstituted in phospholipid nanodiscs. FtsE-His binds an ATP analog with an affinity comparable to that of ATPase subunits of ABC transporters, and FtsE-His preparations have a low, detectable ATPase activity. However, attempts to detect complexes of purified FtsX-His, FtsE-His, and PcsB-His or coexpressed tagged FtsX and FtsE were not successful with the constructs and conditions tested so far. In working with nanodiscs, we found that PcsB-His has an affinity for charged phospholipids, mediated partly by interactions with its coiled-coil domain. Together, these findings represent first steps toward reconstituting the FtsEX:PcsB complex biochemically and provide information that may be relevant to the assembly of the complex on the surface of pneumococcal cells.

Nonfouling NTA-PEG-Based TEM Grid Coatings for Selective Capture of Histidine-Tagged Protein Targets from Cell Lysates.

We report the preparation and performance of TEM grids bearing stabilized nonfouling lipid monolayer coatings. These films contain NTA capture ligands of controllable areal density at the distal end of a flexible poly(ethylene glycol) 2000 (PEG2000) spacer to avoid preferred orientation of surface-bound histidine-tagged (His-tag) protein targets. Langmuir-Schaefer deposition at 30 mN/m of mixed monolayers containing two novel synthetic lipids-1,2-distearoyl-sn-glycero-3-phosphoethanolamine-N-[(5-amido-1-carboxypentyl)iminodiacetic acid]polyethylene glycolamide 2000) (NTA-PEG2000-DSPE) and 1,2-(tricosa-10',12'-diynoyl)-sn-glycero-3-phosphoethanolamine-N-(methoxypolyethylene glycolamide 350) (mPEG350-DTPE)-in 1:99 and 5:95 molar ratios prior to treatment with a 5 min, 254 nm light exposure was used for grid fabrication. These conditions were designed to limit nonspecific protein adsorption onto the stabilized lipid coating by favoring the formation of a mPEG350 brush layer below a flexible, mushroom conformation of NTA-PEG2000 at low surface density to enable specific immobilization and random orientation of the protein target on the EM grid. These grids were then used to capture His6-T7 bacteriophage and RplL from cell lysates, as well as purified His8-green fluorescent protein (GFP) and nanodisc solubilized maltose transporter, His6-MalFGK2. Our findings indicate that TEM grid supported, polymerized NTA lipid monolayers are capable of capturing His-tag protein targets in a manner that controls their areal densities, while efficiently blocking nonspecific adsorption and limiting film degradation, even upon prolonged detergent exposure.

In vitro reassembly of the ribose ATP-binding cassette transporter reveals a distinct set of transport complexes.

Bacterial ATP-binding cassette (ABC) importers are primary active transporters that are critical for nutrient uptake. Based on structural and functional studies, ABC importers can be divided into two distinct classes, type I and type II. Type I importers follow a strict alternating access mechanism that is driven by the presence of the substrate. Type II importers accept substrates in a nucleotide-free state, with hydrolysis driving an inward facing conformation. The ribose transporter in Escherichia coli is a tripartite complex consisting of a cytoplasmic ATP-binding cassette protein, RbsA, with fused nucleotide binding domains; a transmembrane domain homodimer, RbsC2; and a periplasmic substrate binding protein, RbsB. To investigate the transport mechanism of the complex RbsABC2, we probed intersubunit interactions by varying the presence of the substrate ribose and the hydrolysis cofactors, ATP/ADP and Mg(2+). We were able to purify a full complex, RbsABC2, in the presence of stable, transition state mimics (ATP, Mg(2+), and VO4); a RbsAC complex in the presence of ADP and Mg(2+); and a heretofore unobserved RbsBC complex in the absence of cofactors. The presence of excess ribose also destabilized complex formation between RbsB and RbsC. These observations suggest that RbsABC2 shares functional traits with both type I and type II importers, as well as possessing unique features, and employs a distinct mechanism relative to other ABC transporters.

The increasingly complex mechanism of HMG-CoA reductase.

HMG-CoA reductase (HMGR) is the target of statins, cholesterol-lowering drugs prescribed to millions of patients worldwide. More recent research indicates that HMGR could be a useful target in the development of antimicrobial agents. Over the last seven decades, researchers have proposed a series of increasingly complex reaction mechanisms for this biomedically important enzyme. The maturation of the mechanistic proposals for HMGR have paralleled advances in a diverse set of research areas, such as molecular biology and computational chemistry. Thus, the development of the HMGR mechanism provides a useful case study for following the advances in state-of-the-art methods in enzyme mechanism research. Similarly, the questions raised by these mechanism proposals reflect the limitations of the methods used to develop them. The mechanism of HMGR, a four-electron oxidoreductase, is unique and far more complex than originally thought. The reaction contains multiple chemical steps, coupled to large-scale domain motions of the homodimeric enzyme. The first proposals for the HMGR mechanism were based on kinetic and labeling experiments, drawing analogies to the mechanism of known dehydrogenases. Advances in molecular biology and bioinformatics enabled researchers to use site-directed mutagenesis experiments and protein sequencing to identify catalytically important glutamate, aspartate, and histidine residues. These studies, in turn, have generated new and more complicated mechanistic proposals. With the development of protein crystallography, researchers solved HMGR crystal structures to reveal an unexpected lysine residue at the center of the active site. The many crystal structures of HMGR led to increasingly complex mechanistic proposals, but the inherent limitations of the protein crystallography left a number of questions unresolved. For example, the protonation state of the glutamate residue within the active site cannot be clearly determined from the crystal structure. The differing protonation state of this residue leads to different proposed mechanisms for the enzyme. As computational analysis of large biomolecules has become more feasible, the application of methods such as hybrid quantum mechanics/molecular mechanics (QM/MM) calculations to the HMGR mechanism have led to the most detailed mechanistic proposal yet. As these methodologies continue to improve, they prove to be very powerful for the study of enzyme mechanisms in conjunction with protein crystallography. Nevertheless, even the most current mechanistic proposal for HMGR remains incomplete due to limitations of the current computational methodologies. Thus, HMGR serves as a model for how the combination of increasingly sophisticated experimental and computational methods can elucidate very complex enzyme mechanisms.

A novel role for coenzyme A during hydride transfer in 3-hydroxy-3-methylglutaryl-coenzyme A reductase.

In this study, we take advantage of the ability of HMG-CoA reductase (HMGR) from Pseudomonas mevalonii to remain active while in its crystallized form to study the changing interactions between the ligands and protein as the first reaction intermediate is created. HMG-CoA reductase catalyzes one of the few double oxidation-reduction reactions in intermediary metabolism that take place in a single active site. Our laboratory has undertaken an exploration of this reaction space using structures of HMG-CoA reductase complexed with various substrate, nucleotide, product, and inhibitor combinations. With a focus in this publication on the first hydride transfer, our structures follow this reduction reaction as the enzyme converts the HMG-CoA thioester from a flat sp(2)-like geometry to a pyramidal thiohemiacetal configuration consistent with a transition to an sp(3) orbital. This change in the geometry propagates through the coenzyme A (CoA) ligand whose first amide bond is rotated 180° where it anchors a web of hydrogen bonds that weave together the nucleotide, the reaction intermediate, the enzyme, and the catalytic residues. This creates a stable intermediate structure prepared for nucleotide exchange and the second reduction reaction within the HMG-CoA reductase active site. Identification of this reaction intermediate provides a template for the development of an inhibitor that would act as an antibiotic effective against the HMG-CoA reductase of methicillin-resistant Staphylococcus aureus.

Crystallographic analysis of the ENTH domain from yeast epsin Ent2 that induces a cell division phenotype.

Epsins are eukaryotic, endocytic adaptor proteins primarily involved in the early steps of clathrin mediated endocytosis. Two epsins exist in Saccharomyces cerevisiae, Ent1 and Ent2, with single epsin knockouts being viable, while the double knockout is not. These proteins contain a highly conserved Epsin N-terminal homology (ENTH) domain that is essential for cell viability. In addition, overexpression of the ENTH domain of Ent2 (ENTH2) was shown to play a role in cell division by interacting with the septin organizing, Cdc42 GTPase activating protein, Bem3, leading to increased cytokinesis failure. In contrast, overexpression of the ENTH domain of Ent1 (ENTH1) does not affect cytokinesis, despite being 75% identical to ENTH2. An ENTH2(N112D, S114E, E118Q) mutant that switches residues in loop 7 to those found correspondingly in ENTH1 was incapable of inducing the cytokinesis phenotype. In order to better understand the role of loop 7 in the ENTH2-induced phenotype at a molecular level, X-ray crystallography was used to elucidate the structures of yeast ENTH2(WT) and ENTH2(DEQ). Our results indicate that mutations did not affect the conformation of loop 7, but rather introduce an increased negative charge on a potential interaction interface. Morphological analysis of cells overexpressing ENTH2 loop 7 mutants showed that the cytokinesis failure phenotype was abolished by the single mutants N112D, E118Q, and to a lesser extent by S114E. Taken together, our results indicate that the interaction surface that contains loop 7 and the specific nature of these residues are crucial for ENTH2 involvement in cytokinesis. This research provides insight into a molecular mechanism by which ENTH2, but not ENTH1, overexpression in yeast leads to cell division defects. Structural data of WT and mutant ENTH2 domains along with in vivo phenotypic analysis of ENTH2 overexpressing cells indicate that the biochemical nature of three loop 7 residues is crucial for its role in cytokinesis.

Molecular modeling of the reaction pathway and hydride transfer reactions of HMG-CoA reductase.

HMG-CoA reductase catalyzes the four-electron reduction of HMG-CoA to mevalonate and is an enzyme of considerable biomedical relevance because of the impact of its statin inhibitors on public health. Although the reaction has been studied extensively using X-ray crystallography, there are surprisingly no computational studies that test the mechanistic hypotheses suggested for this complex reaction. Theozyme and quantum mechanical (QM)/molecular mechanical (MM) calculations up to the B3LYP/6-31g(d,p)//B3LYP/6-311++g(2d,2p) level of theory were employed to generate an atomistic description of the enzymatic reaction process and its energy profile. The models generated here predict that the catalytically important Glu83 is protonated prior to hydride transfer and that it acts as the general acid or base in the reaction. With Glu83 protonated, the activation energies calculated for the sequential hydride transfer reactions, 21.8 and 19.3 kcal/mol, are in qualitative agreement with the experimentally determined rate constant for the entire reaction (1 s(-1) to 1 min(-1)). When Glu83 is not protonated, the first hydride transfer reaction is predicted to be disfavored by >20 kcal/mol, and the activation energy is predicted to be higher by >10 kcal/mol. While not involved in the reaction as an acid or base, Lys267 is critical for stabilization of the transition state in forming an oxyanion hole with the protonated Glu83. Molecular dynamics simulations and MM/Poisson-Boltzmann surface area free energy calculations predict that the enzyme active site stabilizes the hemithioacetal intermediate better than the aldehyde intermediate. This suggests a mechanism in which cofactor exchange occurs before the breakdown of the hemithioacetal. Slowing the conversion to aldehyde would provide the enzyme with a mechanism to protect it from solvent and explain why the free aldehyde is not observed experimentally. Our results support the hypothesis that the pK(a) of an active site acidic group is modulated by the redox state of the cofactor. The oxidized cofactor and deprotonated Glu83 are closer in space after hydride transfer, indicating that indeed the cofactor may influence the pK(a) of Glu83 through an electrostatic interaction. The enzyme is able to catalyze the transfer of a hydride to the structurally and electronically distinct substrates by maintaining the general shape of the active site and adjusting the electrostatic environment through acid-base chemistry. Our results are in good agreement with the well-studied hydride transfer reactions catalyzed by liver alcohol dehydrogenase in calculated energy profile and reaction geometries despite different mechanistic functionalities.

Specificity of HCPTP variants toward EphA2 tyrosines by quantitative selected reaction monitoring.

EphA2 receptor tyrosine kinase and the human cytoplasmic protein tyrosine phosphatase (HCPTP) are overexpressed in a number of epithelial cancers. Overexpressed EphA2 in these cancers shows a significant decrease in phosphotyrosine content which results in suppression of receptor signaling and endocytosis and an increase in metastatic potential. The decreased phosphotyrosine content of EphA2 has been associated with decreased contact with its ligand, ephrin A1 and dephosphorylation by HCPTP. Potential specificity of the two HCPTP variants for tyrosines on EphA2 has not been investigated. We have used a mass spectrometry assay to measure relative rates of dephosphorylation for the two HCPTP variants at phosphotyrosine sites associated with control of the EphA2 kinase activity or interaction with downstream targets. Our results suggest that although both variants dephosphorylate the EphA2 receptor, the rate and specificity of dephosphorylation for specific tyrosines are different for HCPTP-A and HCPTP-B. The SAM domain tyrosine Y960 which has been implicated in downstream PI3K signaling is dephosphorylated exclusively by HCPTP-B. The activation loop tyrosine (Y772) which directly controls kinase activity is dephosphorylated about six times faster by HCPTP-A. In contrast, the juxtamembrane tyrosines (Y575, Y588 and Y594) which are implicated in both control of kinase activity and downstream signaling are dephosphorylated by both variants with similar rates. This difference in preference for dephosphorylation sites on EphA2 not only illuminates the different roles of the two variants of the phosphatase in EphA2 signaling, but also explains why both HCPTP variants are highly conserved in most mammals.

Use of selected reaction monitoring data for label-free quantification of protein modification stoichiometry.

Quantification of protein and PTM abundance in biological samples is an important component of proteomic studies. Label-free methods for quantification using MS are attractive because they are simple to implement and applicable to any experimental system. We demonstrate that PTM stoichiometry can be accurately measured using label-free quantification and selected reaction monitoring. Use of selected reaction monitoring is advantageous with complex biological samples and we show this approach can be used to quantify multiple PTMs independently on a single peptide.